Nonequilibrium fluctuations in a resistor 
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In small systems where relevant energies are comparable to thermal agitation, fluctuations are 
of the order of average values. In systems in thermodynamical equilibrium, the variance of these 
fluctuations can be related to the dissipation constant in the system, exploiting the Fluctuation- 
Dissipation Theorem (FDT). In non-equilibrium steady systems. Fluctuations Theorems (FT) addi- 
tionally describe symmetry properties of the probability density functions (PDFs) of the fluctuations 
of injected and dissipated energies. We experimentally probe a model system: an electrical dipole 
driven out of equilibrium by a small constant current 7, and show that FT are experimentally ac- 
cessible and valid. Furthermore, we stress that FT can be used to measure the dissipated power 
V — RI^ in the system by just studying the PDFs symmetries. 

PACS numbers: 05.40.-a, 05.70.-a, 07.50.-e, 84.30. Bv 



Introduction 

In the last decade, Fluctuation Theorems (FT) P, Q 
appeared in nonequilibrium statistical physics. These 
new theorems relate the asymmetry of fluctuations of 
energies (or powers) with the dissipated power required 
to maintain a nonequilibrium steady state of a system. 
Thus FT give a measure of the distance from equilibrium. 
Indeed, FT can be extremely useful to probe nonequilib- 
rium states in nanophysics and biophysics where energies 
involved are typicaly small and thermal agitation cannot 
be neglected: in those systems, the variance of fluctua- 
tions is of the order of the average values. Moreover, stan- 
dard Fluctuation-Dissipation Theorem (FDT) 0,1111 is 
derived for equilibrium systems only and so it may fail to 
describe states far from equilibrium. FT are a generaliza- 
tion of FDT out of — possibly far from — equilibrium. 
In order to safely apply FT in complex systems — e.g. in 
biophysics or nanotechnologies — it is important to test 
them in all simple cases where theoretical predictions can 
be accurately checked. In spite of the large theoretical 
effort, just a few experiments have been conducted on 
this topic 0, 0|- For this reason, we test in this letter 
their use on a simple electrical system. 

Electrical systems are particuliarly interesting, because 
they are one of the first where FTD was formalized: 
Johnson H and Nyquist related equilibrium fluctu- 
ations of voltage U across a dipole with the resistive part 
of this dipole. Moreover, all parameters in the setup can 
easily be controled. Thus, out of equilibrium, electrical 
circuits are a key example to stress the differences be- 
tween FT and FDT, and to verify the vahdity of FT. 

Our system is an electrical dipole constituted of a resis- 
tance R in parallel with a capacitor C (Fig.^. We drive 
it out of equilibrium by making a constant current / flow 
in it. Noting fcs the Boltzmann constant and T the abso- 
lute temperature, the injected power is typically of some 
k^T per second, of the same order as in biophysics or 
nanoscale physics experiments. This represent the fun- 
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FIG. 1: Model circuit : an electrical dipole is composed of 
a resistive part R and a capacitive part C. Due to thermal 
fluctuations of charges positions, a fluctuating voltage U is 
observed. We drive the system away from equilibrium by 
imposing a constant flux of electrons, via a constant current 
I. 



damental case of a system in contact with two different 
(electrons-) reservoirs, one of the simplest and most fun- 
damental problems of nonequilibrium physics . U sing 
a powerful analogy with a forced Langevin equation , 
a precise formulation of FT was recently given [12] in this 
case. Nyquist FDT as well as FT are checked experimen- 
tally, by looking at the injected power and the dissipated 
heat in the system. 
Experimental setup 

The circuit we use is composed of a resistor in parallel 
with a capacitor, as depicted on Fig.l. The resistance 
is a standard metallic one of nominal value R = 9.52 
Mf2. In parallel, we have an equivalent capacitor of value 
C = 280 pF. This accounts for the capacitance of the all 
set of coaxial connectors and cables that we used. The 
time constant of the circuit is tq = RC = 2.67ms. Using 
a 50 GQ resistance, we inject in the circuit a constant 
current / ranging from to 6 x 10^^'^ A. This current 
corresponds to an injected power I^R ranging from to 
1000/cBr/s. Experiments are conducted at room temper- 
ature T = 300 K. The typical values of the injected en- 
ergy for, e.g., / — 1.4 X 10~^^ A and r = IOtq are of order 



2 



of a few hundreds of k^T, which is small enough to en- 
sure that the resistance is not heating; expected changes 
of temperature are estimated to be less than 10^^* K 
over a one hour experiment. Moreover, the resistance is 
thermostated: all the heat dissipated by Joule effect is 
absorbed by the thermal bath. 

The fluctuating voltage U across the dipole is measured 
with a resolution of 10~^^ V sampled at 819.2 Hz. This 
is achieved by first amplifying the signal by 10'*, using a 
FET amplifier, with a 4 GO input impedance, a voltage 
noise level of 5 nV/VHz and a current noise of 10~^^ 
A/^/Hz. The signal is then digitized with a 24-bits data 
acquisition card at frequency 8192 Hz, and decimated by 
averaging ten consecutive points. 

Fluctuation-Dissipation theorem 

The electrical dipole of Fig. ^ is a pure resistance R 
in parallel with a capacitor C; its complex impedance 
reads Z{f) = 1/ [1/ R + i2nRC f) where = -1 and / 
is the frequency. The effective dissipative part is the real 
part 3?(Z) of Z . It was experimentally observed by John- 
son 1^, and then demonstrated by Nyquist that the 
potential difference U across the dipole fluctuates with a 
stationary power spectral density S{f)df such that 



Sim = ^knm{z)df 



(1) 



At equilibrium / = 0, in average no current is flowing 
in the circuit, and mean U is zero. Integrating over all 
positive frequencies, one gets the variance of U: 



knTR knT 



To 



c 



(2) 



Eqs. (P) and © are the expressions of the Fluctuation- 
Dissipation theorem (FDT) for electrical circuits. 

The exact value of our capacitance was determined by 
fitting the power density spectrum of equilibrium fluctu- 
ations (at imposed / = A) by a Lorentzian low-pass 
transfer function (eq.(E3)), as illustrated on Fig.|2l3. Ap- 
plication of FDT leads with a very good accuracy to the 
determination of i?, in perfect accordance with the mea- 
sured nominal value (Fig. 2). When / = 1.4 x 10~^^ A, 
we found the same power spectral density for U , and per- 
forming the same treatments gave the same estimates of 
R and C. We therefore conclude that FDT is still holding 
in our system driven out of equilibrium. 

Fluctuation theorems 

The power injected in the circuit is Vm = UI, but 
only the resistive part dissipates, so the dissipated power 
is "Pdiss = f^*R, where zr is the current flowing in the 
resistor (Fig. P). As already noted 0|, in average, one 
expects (■Pin) = ('Pdiss) = where the brackets stand 
for time average over sufficiently long times compared 
to tq. This is very well checked in our experiment. Vm 
and Pdiss fluctuate in time because U itself is fluctuating. 
If one assumes that fluctuations of U have a Gaussian 
distribution, which is the case at equilibrium when / = 0, 
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U,(t) = ^ ff f^^ U(t') dt' 
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FIG. 2: (a): For / = A, Johnson-Nyquist noise has a 
Gaussian distribution. Relation ||5J is verified, (b); The noise 
is white up to the cutoff frequency fo of the RC dipole. We use 
FDT to extract from the energy of the noise the value of the 
resistive, dissipative part of the circuit R: for low frequency, 
the noise level spectral density is constant, equal to y/AksTR 
in a band [/; / -I- d/]. A Lorentzian fit of the spectrum (thin 
line) additionally gives to — RC. 



then Vin has also a Gaussian distribution, because / is 
constant. On the contrary, in fluctuates, as we see from 
Kirchoff's laws: 

I = lR + C'^, SO T'in - T'diss + , (3) 

dt 2 dt 

and therefore, the probability distribution of T'diss is 
not Gaussian [l2|. It is worth noting that for large cur- 
rent /, some orders of magnitude larger than the one we 
use. Via and Pdiss will be much larger than the conser- 
vative part ^ ^gj- and therefore the probability distribu- 
tions of both the injected and dissipated power will be 
practically equal and usually Gaussian, as it is expected 
in macroscopic systems. 

We can {g)r{t) 
value of a function g over a time r 



i /j*^^ g{t')dt' the time-averaged 
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FIG. 3: Histograms of Wt and Qt when a current 7 = 1.4 x 
10~^^A is flowing in the dipole. This corresponds to {Wt) = 
{Qt) = Tp with V = RI'^ ^ 45fcBT/s. 

Reasoning with energies instead of powers, we define 
Wr{t) = T{Vin)T{t), the energy injected in the circuit 
during time r, analogous to the work performed on the 
system (positive when received by the system). In the 
same way, we write Qrit) = T(7^diss)r (0; the energy dis- 
sipated by Joule effect during time r, analogous to the 
heat dissipated by the system (positive when given by the 
system to the thermostat). We used values of r spanning 
from fractions of tq up to hundreds of tq. 

For a given value of /, we measure U (t) and compute 
Vin and Pdiss- We then build the probability density 
functions of cumulated variables Wr and Qr using 10^ 
points; their typical distributions are plotted on Fig. |31 
As expected, fluctuations of Wt are Gaussian for any r 
whereas heat fluctuations are not for small values of t. 
They are exponential for small r: large fluctuations of 
heat Qt are more likely to occur than large fluctuations 
of work Wt- 

Then we study the normalized symmetry function: 

c t ^ - '^bT pjET = g) 

Mr,a) = -^ln^^^— (4) 



FIG. 4: Normalized symmetry functions Sw and Sq for Wt 
and Qt when I — 1.4 x lO'^^'^A (same conditions as Fig. |^. 
(a): For any r, Sw is a linear function of a = Wt/{Wt)- 
Straight lines are theoretical predictions from [l^ . For r 
CO, Sw{a) tends to have a slope 1. (b): On the contrary, 
Sq (a) tends to a limit function (continuous black curve) which 
is a straight line of slope 1 for a < 1 only. 

where Et stands for either Wt/{Wt) or Qt/{Qt)- If the 
Fluctuation Theorem for the work Wt holds, then one 
should have 0, 0| , for large enough r, the relationship 
IuHt ^ oc fw{T, a) = a. In contrast, if the Fluctuation 
Theorem for the heat Qt holds, then for t oo, the 
asymptotic values of SqIt^q) are 5*0^(0) = a for a < 1, 
SQ'{a) = 2 for a > 3, and there is a continuous parabolic 
connection for 1 < a < 3 that has a continuous deriva- 
tive [lll[]3 |. From histograms of Fig.O we compute the 
symmetry functions Sw{T,a) and SQ{T,a) (Fig.0J. 

Work fluctuations First, for any given t we checked 
that the symmetry function Sw {t, a) is linear in a 
(Fig. EJi). We measured the corresponding proportion- 
ality coefficient aw{T) such that SwiT^a) — (7w{t)o-- 
This coefficient crw(r) tends to 1 when r is increased 
(see Fig. El . 

Heat fl,uctuations We found that Sq (r, a) is linear in 
a only for a < 1, as expected 0,^2- Again, as r ^ 00, 
the limit slope of the symmetry function is 1 whereas for 
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FIG. 5: (a): dependence on r of the slope of Sw in Fig. 2] 
Slope is converging from above, (b): distance between Sq{t) 
for finite t and theoretical prediction of Sq for infinite r, for 
several values of a = Q/{Q)- for a < 3, convergence is from 
below, whereas it is from above for a > 3. 



a > 3, Sq{t, a) tends to two. 

Asymptotic symmetry functions and convergence 
In [n, 12], expressions for the convergence towards the 
asymptotic Umits S^{a) and SQ{a) are given in terms 
of T. We can check these predictions with our data. The 
convergence for the work Wr is very well reproduced by 
these predictions, as can be seen on Fig. continuous 
straight lines are theoretical predictions for small values 
of T, using eqs. (9) and (10) from with no adjustable 
parameters. On Fig. |5ji the slope criy(r) of the experi- 
mental symmetry function is plotted. The continuous 
line is the prediction of jTj |. which perfectly agrees with 
our data. 

For the heat, we distinguish two regimes for the con- 
vergence towards the asymptotic symmetry function. For 
a = Qr/ (Q) < 3, we find that when r is increased, sym- 
metry functions are converging to the asymptotic func- 
tion from below (Fig. 0)1, which is the opposite of what 
is observed for the work. On the contrary, for a > 3, 
convergence to the asymptotic function is from above, 
thus enhancing the peculiarity of the point a = 3. On 
Fig. we have plotted the evolution of Sq (t, a) versus 
T for several fixed values of a. The convergence from 
above for a > 3 and from below for a < 3 is clear. For 
increasing values of r, only smaller and smaller values of 
a are accessible because of the averaging process. There- 
fore the accessible values of a are quickly lower than 1, 



and only the linear part of Sq (r, a) can be experimen- 
tally tested. Nevertheless, for intermediate time scales r 
we see in Fig. ^ that the data converge towards the theo- 
retical asymptotic nonlinear symmetry function [ill 
(smooth curve in Fig. ^). 

We observed that the convergence reproduced on FigJE] 
depends on the injected current /, as pointed out in |l2l |. 
Other experiments with a larger current {P = ISG/cbT/s) 
give a faster convergence; corresponding results will be 
reported elsewhere. 

Conclusions 

We have shown experimentally that the asymmetry of 
the probability distribution functions of work and heat in 
a simple electrical circuit driven out of equilibrium by a 
fixed current /, is linked to the averaged dissipated power 
in the system. The recently proposed Fluctuation The- 
orems for first order Langevin systems are then experi- 
mentally confirmed. Exploiting formula Q), FT can be 
used to measure an unknown averaged dissipated power 
V ~ limi-^oo — by using only the symmetries of the 
fluctuations, i.e. computing Sw or Sq and measuring 
their asymptotic slope. 

We operated with energies of order of in order 
to have strong fluctuations compared to the averaged 
values. It is worth noting that as the driving current 
is increased up to macroscopic values, the fluctuations 
become more and more negligible, therefore Fluctuation 
Theorems become harder and harder to use, and there- 
fore less relevant. 
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